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Light scalar fields with only gravitational strength couplings are typically present in UV complete 
theories of physics beyond the Standard Model. In the early universe it is natural for these fields 
to dominate the energy density, and their subsequent decay - if prior to BBN - will typically yield 
some dark matter particles in their decay products. In this paper we make the observation that a 
Non-thermal WIMP 'Miracle' may result: that is, in the simplest solution to the cosmological moduli 
problem, non-thermally produced WIMPs can naturally account for the observed dark matter relic 
density. Such a solution may be generic in string theory compactifications. 

PACS numbers: 



I. INTRODUCTION 

For several decades, a compelling theoretical picture 
for dark matter has been developed and widely applied, 
based on what is now called the Thermal WIMP 'Mira- 
cle'. In this picture, the early Universe is very hot and 
dense and essentially all particle species are in thermal 
and chemical equilibrium. As the Universe expands and 
cools to a temperature of order the mass of the dark 
matter particle Xj the annihilations of x cease to be effi- 
cient at reducing the number of particles compared to the 
cosmic expansion and a 'freeze-out' occurs. The result- 
ing relic number density of x-particles then depends only 
on the ratio of the annihilation cross-section of x ^-nd 
the Hubble scale near the freeze-out temperature. For 
Weakly Interacting Massive Particles (WIMPs), i.e. par- 
ticles with electroweak interactions and masses, the 'mir- 
acle' is that this ratio is in good agreement with that de- 
duced from astrophysical and cosmological experiments 
(see [l| for reviews). 

Whilst this is an extremely compelling idea, it can of- 
ten be difficult to implement in practice, with the relic 
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density predicted by theory disagreeing with the data by 
a couple orders of magnitude in both directions. That is, 
the thermal WIMP 'miracle' faces significant challenges 
when directly confronted with precision data. It is not 
clear at present if this should be viewed as a failure of 
the theoretical models or as a phenomenological guide to 
select particular classes of models. For example, in su- 
persymmetric models, where the thermal relic idea has 
perhaps been most extensively explored, most of the pa- 
rameter space yields an incorrect dark matter density, 
prompting attempts to look at special regions [5| . 

Furthermore, by considering the high energy behavior 
(UV completion) of phenomenologically based models, 
it becomes less straightforward to motivate cosmologies 
in which the Universe is in thermal equilibrium prior to 
Big Bang Nucleosynthesis (BBN). This is partly because 
many UV completions lead to the inclusion of moduli - 
neutral scalar fields which couple to matter only gravita- 
tional - and for a wide range of masses these moduli will 
typically evolve to dominate the cosmic energy density 
prior to BBN. 

As a result, the decay of the moduli must give rise 
to temperatures above that of Big Bang Nucleosynthesis 
(BBN) - a few MeV - or the predictions of BBN will be ru- 
ined. This is the cosmological moduli problem (formerly 
known as the Polonyi problem) 0, 0, H, 0] and typically 
requires that the moduli masses are at least 10 TeV or 
greater. Of course, exceptions are possible. A period 
of low-scale inflation could dilute the moduli density 0] ■ 
Another possibility is if signigicant energy is not stored 
in the fields. This could be accomplished by dynamics. 
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or if the moduli are stabilized initially near points of en- 
hanced symmetry [1, 0. However, in this case, the need 
for a perturbative low energy theory necessarily requires 
at least one modulus to remain unprotected, and thus 
a substantial contribution to the energy density should 
still be expected Q . We will not discuss these interesting 
possibilities further here. 

Moduli which evade the cosmological moduli problem 
have masses which are typically greater than that of 
WIMP dark matter candidates, so their decays will nec- 
essarily result in the production of some dark matter. In 
fact, the yield of dark matter particles can be quite large, 
since moduli tend to couple to the matter sector univer- 
sally with gravitational strength couplings (~ l/rrip) and 
the branching ratios to stable particles are expected to be 
substantial. Thus, the density of dark matter at the time 
of moduli decay will typically be a sizable fraction of the 
energy density of the moduli themselves. Non-thermally 
produced dark matter has been considered previously by 
many authors [13, [S Q [13 ■ 

In this paper we will attempt to elaborate the con- 
nection between the moduli problem and non-thermally 
produced dark matter. In the absence of any strong guid- 
ance as to how one might proceed from fundamental the- 
ory, we will take the viewpoint of a generic supergravity 
theory containing moduli. That is, we will simply take 
an effective field theory approach to the general problem 
of moduli coupled to matter in a supersymmetric frame- 
work. A priori though, in order to proceed, one still 
requires a clue about the mass spectrum of superpart- 
ners, moduli and the dark matter candidate. In fact, the 
effective theory itself provides such a clue. 

To see this, note that in a generic supergravity theory 
the moduli masses are of order 7713/2 (the gravitino mass) 
and, furthermore, all scalars receive mass term contribu- 
tions of the same order This includes all squarks 
and sleptons. Thus, since m3/2 is required to be at least 
10 TeV because of the moduli problem, squarks and slep- 
tons are at least as massive. Further, as we will argue 
below, gauginos need not be as massive as the squarks 
and sleptons in such a generic supergravity theory. On 
the other hand, in the absence of special symmetries, hig- 
gsinos also generically have a mass of order TO3/2 by the 
Giudice-Masiero mechanism. If we thus consider the sim- 
plest model within this framework - the minimal super- 
symmetric standard model (MSSM) coupled to moduli, 
these arguments suggest that x is a (neutral) gaugino, 
i.e. a Bino, Wino or a mixture thereof. 

Rather surprisingly, within the class of models with 
spectra suggested by the above arguments, there exists a 
Non-thermal WIMP 'Miracle' under very general condi- 
tions. In other words, by solving the cosmological moduli 
problem, one automatically obtains a consistent solution 
to the basic dark matter problem with WIMPs when the 
moduli masses are lO's of TeV.. However, as will be seen, 
compared to the thermal case, the WIMP annihilation 
cross-section is larger, thereby suggesting regions of pa- 
rameter space previously avoided based on the thermal 



dark matter picture. 

Examples of models with the sorts of spectra consid- 
ered here have been considered previously in the context 
of AMSBjlJ] and more recently in the context of M the- 
ory [HI . Il2j| where the non-thermal production of dark 
matter was emphasized (see also [ISJ). These examples 
show that with additional, well motivated assumptions, 
that the arguments of this paper can be sharpened and 
detailed models can be constructed. We emphasize that 
many of the ideas given here were inspired by the seminal 
work of [Hi . 



II. MODULI MASSES AND SUPERSYMMETRY 
BREAKING 

The moduli fields must enter the supergravity poten- 
tial. If not, they would remain massless and not gain 
the vevs required to explain, e.g. the value of the fine- 
structure constant aem or the Yukawa couplings in the 
Standard Model. If V is the potential for all scalar fields, 
and the only supersymmetry breaking physics is gener- 
ated at the scale of order VF, then the moduli masses 
will generically also be of order the gravitino mass as we 
now review. 

The supergravity potential evaluated in the vacuum is 
of the form 



(1) 



where Fi are the vacuum values of the SUSY-breaking 
.F-terms which are to be summed over all scalars (f)i with 
non-trivial _F-terms, K is the Kahler potential, and 7713/2 
is the gravitino mass. Hence, because the observed vac- 
uum energy today is so small. 



^3/2 



F/nip 



(2) 



where F is of order the dominant _F-term. This must 
be true, regardless of the mechanism of supersymmetry 
breaking. In fact, F/nip sets the typical mass scale for 
all scalar fields appearing in y in a generic supergravity 
theory. For instance, since 

dW dK 

where W is the superpotential, there are terms in V of 
order 



V ~ KiK' 



(4) 



If we consider the terms in K (or in KiK^ ) of order 0^0* 
we obtain contributions to V like 



V . 



^^<l>i — 2" 



(5) 



^3/2"'p 



In the vacuum \W\'^ 

mass terms for of order TO3/2. Hence, all scalars re- 
ceive contributions to their masses of order 7773/2 in a 
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generic supergravity theory. This includes moduh as well 
as matter scalars such as higgses, squarks and sleptons ^. 

Exceptions exist to the above general statements. In 
the language of low energy supersymmetry, these would 
correspond to very special vacuum Kahler and superpo- 
tentials for the moduli and matter fields. For instance, if 
the moduli dynamics is i?-symmetry preserving, then the 
moduli potential will not break supersymmetry and the 
moduli could, in principle, obtain large masses and vevs 
independent of the value of the gravitino mass, whose 
value is set by additional i?-breaking dynamics at an- 
other scale. We will not consider such exceptions further, 
except in the conclusions, and will adopt the viewpoint 
of the 'generic supergravity Lagrangian'. 

Notice also that, if there were a bare, large mass term 
in W for the moduli, the fact that the moduli get vevs 
will give W a vev in general and therefore contributes to 
7713/2 thereby connecting again the moduli and gravitino 
mass. Again, to avoid this requires i?-symmetric moduli 
dynamics. 

A. Gauge Mediation 

In gauge mediation the dominant F-term \/F < 10^° 
GeV, though typically F is taken to be much smaller 
than what the upper limit (of high scale gauge mediation) 
suggests. Hence, the gravitino mass and moduli masses 
are of order 

7773/2 ^ ^ F/wLp < lOGeV (6) 

Since theories of gauge mediated supersymmetry 
breaking have \/F between a TeV and 10^° GeV they 
lead to moduli masses in the wide range: 

lO^^eV < 7770 < 10 GeV (7) 

Of course, in gauge mediation, charged scalars get 
much larger corrections to their masses from their in- 
teractions with the 'messengers of SUSY breaking', but 
the moduli do not. 

This range of moduli masses leads to cosmological 
problems. First of all, since the moduli have Planck scale 
suppressed couphngs to matter, this range of masses gives 
moduli whose lifetimes are between 100 and lO''^ years. 
After inflation, when the Hubble scale becomes of order 
7710 the moduli begin to oscillate and will quickly dom- 
inate the energy density of the universe over radiation. 
This range of masses corresponds to temperatures 

TeV < Tosc < 10^° GeV (8) 

which means that such moduli would begin to dominate 
before BEN. Therefore, such models lead to Universes 



^ See [Tsll for examples. 



which are dominated by moduli for very long periods. 
Further, after the moduli decay the Universe is not re- 
heated enough to start BBN with the correct conditions. 
Hence, gauge mediation models coupled to moduli re- 
quire quite special Kahler and super potentials which 
would allow for the moduli to be very massive compared 
to the gravitino. 



B. Moduli Masses in Gravity Mediation 

In gravity mediation, \/F is much higher, e.g. lO^^GeV 
to lO^^GeV and hence 7713/2 is usually taken to be of order 
TeV, since supersymmetry is assumed to solve the hier- 
archy problem. Then, through the typical supergravity 
couplings discussed above, the moduli will also end up 
with TeV scale masses. Such moduli will decay during 
BBN and will typically ruin its successful predictions for 
light element abundances. The gravitino also leads to 
similar problems. This is just the usual moduli problem. 

In gravity mediation, there is a simple solution to these 
problems. One can simply raise the scale of SUSY break- 
ing (Vf) by a factor of a few which raises both 7773/2 and 
7770 by one order of magnitude. This decreases the mod- 
uli lifetime by three orders of magnitude and is consistent 
with BBN occurring just after the moduli have decayed. 
A detailed model in this case has been described in [13 |. 
Notice that the gravitino mass has also been raised above 
the TeV scale, since F is increased. 

Since, as we have discussed above, all scalars, including 
squarks, sleptons and Higgses, will have masses of order 
7773/2 the fine tuning problem between 7773/2 and mz is 
naively much worse in theories with moduli than in the 
usual little Hierarchy problem. 

A crucial point, however, is that the gauginos do not 
have to be as massive as the squarks and sleptons: they 
can easily be lighter if the field whose F-term dominates 
SUSY breaking in the hidden sector is not the field whose 
vev generates the gauge couplings a. In fact, in the 
generic supergravity point of view, there is no reason 
why the gauge coupling function should be dominated 
by the field with the dominant _F-term. This can pre- 
sumably also be understood as the consequence of an 
approximate i?-symmetry in this sector. Here we em- 
phasize that it is quite generic. However, in general this 
approximate i?-symmetry only suppresses the gaugino 
masses but not that of the higgsinos. One could also 
have special i?-symmetries in which both gaugino and 
higgsino masses are suppressed, as in split supersymme- 
try 19']. Hence, we expect that the gauginos are sig- 
nificantly lighter than the squarks and sleptons. In the 
absence of special i?-symmetries, we expect the higgsinos 
to get masses through the Giudice-Masiero mechanism of 
order 7773/2. 

What we learn is - by considering the simplest solution 
to the cosmological moduli problem, we obtain a rough 
picture of the spectrum of BSM particles. Existing mod- 
els which have this kind of spectrum include Anomaly 
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Mediated Supersymmetry Breaking Models (AMSB) [Ij 
and the G2-MSSM [ll. 



C. The Non-thermal WIMP 'Miracle' 

By taking this effective field theory approach we have 
arrived at a picture which suggests that supersymmetry 
breaking is gravity mediated with moduli and gravitino 
masses of order 10 TeV. The squark, slepton and higgsino 
masses are also of order 10 TeV, whereas the gauginos are 
typically lighter. Cosmologically, these moduli inevitably 
dominate the Universe after inflation and up to BBN. 

We can now ask: with the spectrum roughly fixed by 
these arguments, how much dark matter is produced, 
both thermally, and non-thermally by the decays of the 
moduli fields? In other words, is there a connection be- 
tween the moduli and dark matter problems? 

After inflation, when the Hubble expansion becomes 
comparable to the mass of the moduli {m^) they will 
begin to oscillate in their potential. The resulting energy 
density will dilute like normal pressure-less matter and 
they will quickly come to dominate the energy density. 
The oscillations will begin at a temperature determined 
by the moduli mass 



T 



I ^l/2 



(9) 



The resulting condensate will then decay when the ex- 
pansion rate becomes comparable to the decay rate 



H 



(10) 



Given that the moduli dominate before their decays, the 
resulting energy density is (assuming instantaneous de- 
cay) Pdecay ~ T^^TO^ = m^/m^ and so the expected 
number density of dark matter particles will be 



(11) 



We can compare with the critical number density for an- 
nihilations to occur. 



H 



av av m'^{av)' 



(12) 



where av is the self annihilation cross-section times ve- 
locity of the produced dark matter particles. 



[av 



(13) 



the cosmological moduli problem (~ 10 — 100 TeV), we 
find that the number density of produced particles is 
easily large enough for the x particles to annihilate effi- 
ciently, unless the branching ratio into dark matter parti- 
cles is very small. This is unlikely if dark matter particles 
are gauginos, because the gauge coupling is a modulus 
vev and this has an order one coupling to the gauginos. 
They will thus continue to annihilate until their number 
density becomes of order — . This is the non-thermal 
analogue of 'thermal freeze out'. Hence, the non-thermal 
freeze out number density is 



(14) 



where the temperature of the Universe after the decay is 

3/2 



Trh ~ (r^mp) 



1/2 



1/2 



(15) 



The fact that a non-thermal freeze out actually oc- 
curs is by itself an important statement, making the relic 
density more model independent than if it did not occur. 
Assuming no further entropy production in the Universe 
after this stage, the resulting relic abundance is 



; 0.1 X 



100 GoV^ 



10.75 



1/4 



(av) 



100 TeV 



J 



3/2 



(16) 

where ctq — 3 x 10 ^''cm'^s ^ and we have assumed^ 
the constant of proportionality in pop is around a fac- 
tor of ten. For x particles with order 100 GeV mass, 
a 100 TeV moduli mass scale gives excellent agreement 
with the data. Note that, (fT4|) expresses the number den- 
sity at freeze-out as the ratio of the Hubble parameter to 
a particle physics cross-section, just as in the thermal 
case, except now the Hubble scale is evaluated at the 
moduli reheating temperature. This is the non-thermal 
WIMP 'miracle'. With larger moduli masses, smaller 
cross-sections are also possible. However, this separates 
the Electroweak scale from the scale of supersymmetry 
breaking even further and, for masses beyond 1000 TeV 
there are no good dark matter candidates in the spec- 
trum whose cross-section is small enough, leading to too 
little non-thermal dark matter. Therefore, the preferred 
moduli mass scale is roughly 10 to 100 TeV. 

For 0(100) GeV masses, the x particles freeze out at 
temperatures of order a few MeV, so the annihilation 
cross-section is about two to three orders of magnitude 
larger than the thermal case, where freeze out occurs at 
a few GeV. Also note that larger moduli masses dilute 
the final relic density, so that if the dark matter is to be 
explained this way by gauginos, it requires moduli masses 
just beyond the weak scale. 



Taking typical weak scale values for the dark matter 
mass {r-^ 100 GeV) and cross-section 10~^'*cm'^s~i ~ 
10~^GeV~^), and moduli masses in the range to address 



^ Larger values for the proportionality constant were found from 
explicit calculations in |1^ . and smaller values are possible as 
well (see e.g. \V^ ). 
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III. DISCUSSION 

The consistent solution of both the dark matter and 
moduh problems suggests that there is another impor- 
tant scale in nature, of order 10-100 TeV. The picture 
which emerges here is, in some sense, as compelling as the 
thermal case, but has the additional advantage that it is 
consistent with string theory and other frameworks with 
moduli fields, without having to invoke 'unknown dynam- 
ics' which decouple the moduli. The non-thermal dark 
matter picture requires larger annihilation cross-sections 
than the thermal case, thereby suggesting a class of mod- 
els which would not be considered based upon the ther- 
mal relic density picture. Larger cross-sections are also 
helpful to explain the cosmic ray positron excess observed 
by PAMELA. For instance, the neutral Winos (with an- 
nihilation cross-section av « 2.4 x 10"^'' cm^ • s~^) which 
are the dark matter particles ir i [ill [T^ could be a good 
fit to the PAMELA data [HI, |22|. Non-thermally pro- 
duced dark matter particles will in general have a dif- 
ferent phase space distribution than the thermal case; 
in principle this could lead to observable consequences 
for structure formation and the cosmic microwave back- 
ground. Finally, although entropy production by the de- 
cay of moduli also dilutes any pre-existing baryon asym- 
metry, the problem can be naturally solved if either a suf- 
ficiently large initial baryon asymmetr y (s uch as by the 
Affleck-Dine mechanism) is generated [2J|, or if the de- 



cays of the moduli themselves generate both the baryon 
asymmetry and Dark Matter [2J| . 

We also point out that these arguments tend to disfa- 
vor gauge mediated supersymmetry breaking. If gauge 
mediated supersymmetry breaking were discovered at 
the LHC, it would imply that our supergravity vacuum 
is extremely non-generic, in that it allows large moduli 
masses whilst keeping the gravitino very light. An exam- 
ple might be a theory in which the moduli dynamics is 
_R-symmetric whilst supersymmetry breaking is not. This 
might allow a decoupling of 7713/2 and 777^. However, as 
noted in Q, decoupling these scales is very difficult in 
general. 
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